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ABSTRACT 

The heat conduction equation for the neutral atmosphere 
in the region between 120-1500 km is solved using the concept 
of dynamic diffusion. It is shown that the present discrepancy 
in the phase and the diurnal amplitude of the thermospheric 
temperatures as inferred from incoherent back-scatter and 
satellite drag measurements can be fully reconciled by such an 
approach. The numerical solutions show that the phase and the 
diurnal amplitude of the atmospheric density in the region between 
200-500 km are in accord with the satellite drag results; the 
corresponding neutral temperature, however, follows the diurnal 
pattern inferred from the back-scatter technique. The charac- 
teristics of the density variations above 500 km strongly depend 
upon the choice of the hydrogen and helium models and their 
relative concentration. In the region where hydrogen becomes a 
dominant constituent, the density reaches its maximum value 
in the nighttime and its minimum in the day-time. 

The physical basis for ty'e dynamic diffusion is discussed, 
and it is suggested this should replace the conventional static 



diffusion models in computing temperatures from the atmospheric 
densities obtained from satellite drag measurements. 



INTRODUCTION 


The diurnal characteristics of the thermospheric density 
and temperature and their physical interpretations have been 
the subject of considerable discussion during the past few 
years. The satellite drag results which have been the major 
source of information about the neutral atmosphere over the 
past decade seem to indicate that the atmospheric density and 
temperature above 200 km reach their maximum values around 
1400 hours local time almost universally. (Jacchia and Slowey 
1968, See also Priester el. al (1967) for other references). 

The theoretical considerations which are based on the solution 
of the heat conduction equation with solar ultra volet radiation 
as the main energy input , suggest that at equinox the diurnal 
temperature should have its peak value around 1700 hours . In 
order to explain this phase discrepancy of 3 hours, Harris and 
Priester (1962) introduced the concept of a secondary heat 
source to which they ascribed a corpuscular origin. In spite 
of the success of Harris and Priester in bringing their 
theoretical calculations in accord with the satellite drag 
measurements, the concept of this new heat source has been 
seriously questioned in the literature. A number of alterna- 
tives have been suggested to this heat source which include the 
introduction of dynamical effects, coupling with the mesosphere 
and extending the one dimensional calculation of the heat con- 
duction equation to two or three dimensions. (Freeman 1967, 
Dickinson et al . 1968, Volland et al . 1969). All these papers 


3 


succeed to varying degrees in bringing the theoretical calcu- 
lations in accord with the satellite drag results without the 
necessity of an additional heat source. The controversy now 
has taken a new turn. The recent observations on the thermo- 
spheric temperatures inferred from the incoherent back- 
scatter radars suggest the temperature to have its peak value 
around 1700 hours in agreement with the calculations of Harris 
and Priest er without an additional heat source and in sharp 
disagreement with the satellite drag results. (Nisbet 1968, 
Mahajan 1968, McClure 1968, Carru et al . 1968, Waldteufel and 
McClure 1970) . The disagreement between the satellite drag 
and the incoherent back-scatter extends beyond the time of the 
maximum. The minimum temperature inferred from the back-scatter 
measurements occurs about one hour later than the time of the 
minimum inferred from the drag measurements. The diurnal 
amplitude of the temperature inferred from the two sets of 
measurements are also different. The ratio of maximum to mini- 
mum temperature derived from the back-scatter data is about 1.5. 
The corresponding value from the drag measurement is about 1.3. 

The purpose of this paper is to show that the phase and 
the amplitude of the neutral atmosphere derived from the back- 
scatter and drag measurements do not contradict each other 
if is recognized that the basic parameters measured from the 
two experiments are the temperature and the density, respectively. 
The basic problem lies in computing the temperature from the 
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density and vice versa. It is proposed in this paper that the 
concept oi the dynamic diffusion as opposed to the static dif- 
fusion (Jacchia 1964) usually employed in computing the tempera- 
ture from the density leads to a reasonable reconciliation 
between the two sets of the observational results. 

THE P HYSICA L CONCEPT OF DYNAMIC DIFFUSIO N 

The concept of the dynamic diffusion is based on the 
assumption that the atmosphere can be divided into two regions 
according to the relative importance of diffusion or mixing as 
the controlling process. The turbopause which is usually taken 
to be the height above which diffusive equilibrium prevails in 
the atmosphere is assumed here to be also the level separating 
the mixing and diffusive regions. In reality there is no sharp 
boundary between the two regions and this reference level is 
chosen here for convenience only, without implying a new defi- 
nition for the turbopause. Both the density and temperature 
in the two regions undergo diurnal variations as determined by 
the dynamical and heat transfer processes of that region. How- 
ever, the changes in the two regions are not independent of 
each other because of the coupling through a common boundary. 

The solutions of the continuity, momentum and energy balance 
equations must match at this level. Accordingly, the density 
and temperature at the turbopause are not invariant but undergo 
temporal and diurnal variations. Based on this idealized con- 
cept, the particle density n. of the jth constituent of the 

J 

neutral atmosphere above the turbopause may be expressed in the 
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iol lowing form: 
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where k is the Boltzmann constant g the acceleration due to 
gravity; and n (z Q ,t) and T (z Q ,t) are the particle density 
and temperature at the turbopausc level z . 

Equation (1) essentially summarizes the basic concept of 
dynamic diffusion which may be stated as follows: At any 

given instant, the neutral constituents in the upper atmosphere 
above the turbopause are in diffusive equilibrium even though 
the density and temperature at the lower boundary are not in- 
variant. The variations at the lower boundary which are caused 
by the physical processes lower down tend to create transients 
in the medium. However, because of an extremely short time con- 
stant of the atmosphere, these transients die very quickly and 
the transition from one diffusive equilibrium state to another 
is almost instantaneous. Thus, if the variations at the lower 
boundary are known, the diurnal variations in neutral density 
and temperature in the diffusive region can very simply be cal- 
culated by solving the heat conduction equation in conjunction 
with equation (1) . Because of the variability at the lower 
boundary, it is evident that the temperature and density at any 
given height will not be in the same phase and will depend 
strongly on the phase and amplitude of the diurnal variations 
at the lower boundary. We shall return to this point in the 
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later part of this paper while discussing the numerical results 
obtained from the solutions of the heat conduction equation. 

Since the concept o 1 dynamic diffusion and hence* the main 
conclusion of this paper is based on the assumed variability 
at the lower boundary, it i r > important to investigate* the 
nature and extent of these variations. This is considered in 
the next section. 

THE DIURNAL V ARIATIONS BELOW THE TURBOPAt ) SE 

The discussions of this section >vill have a rather qualitative* 
character and are meant mainly to make plausible the possibility 
of appreciable density variations at the turbopause level . The* 
result, expressed in Eqn . (14), will be a formula giving an 
approximate relationship between the density and the temperature 
at the turbopause. Since we are mainly interested in the density, 
we will have to specify the temperature with respect to its 
diurnal amplitude and phase, and we will obtain in an approximate 
way the corresponding density variation. We will simplify our 
problem by assuming that the atmosphere below the turbopause is 
completely mixed and the main dynamical effects of this region 
are controlled by the temperature variations which arc* specified. 
This seems to be a reasonable assumption for the region between 
70 - 120 km where solar X-rays, Lyman a and the ultraviolet 
radiation in the Schumann continuum range piovide the main heat 
source in the atmosphere. The equation of continuity and 
momentum conservation for this region, neglecting the viscosity 
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and the tidal effects, may be written in the* following 
l'orm 
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whore n is the angular velocity (i 7.27 x 10 rad/sec) of 
the earth, r and v art* respect ively the density and the velocity 
of the gas, and p is the scalar pressure 


p n k T (4) 

For a given temperature distribution, (2) and (3) along 
with (4) can be uniquely solved subject to the boundary con- 
ditions. We assume here that the changes in temperature, density 
and pressure are small compared to their mean values. This 
enables us to adopt the usual perturbation technique for solving 
these equations. Thus, if we assume that all the time varying 
quantities can be expanded in a Fourier series, we may write 
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where cr is the fundamental frequency. The quantities p^ r , 
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T lr , and v Jr art* complex, thus containing a phase term. 

Substituting (5) in (2), (3) and (4), the resulting dil- 
Jerential equations can be solved lor lr and v 1 . for a hori- 
zontally stratified atmosphere in which the importance of 
horizontal pressure gradients can be neglected as compared 
to the vertical, lr and Vj rz are simply given by the following 
expressions (See Appendix) . 
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where A and B are constants of integration and are in general 
complex, K r , r QO and H are defined in the appendix. 

The constants A and B can be evaluated by assuming that at 
the lower boundary (z = o) , v^ rz . = o . At the upper boundary 
(z = z Q ) from where the diffusive equilibrium begins we note that 
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corresponding to the definition of the vertical velocity 


as the velocity of isobaric surfaces. Thus, at z z 
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From these constraints on the vertical velocity. A and D 
can be shown to have the following expressions 
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The expression for density variations at the turbopause in 
which we are primarily interested in this problem is obtained by 
substituting z z q in equation (6). From (6), (11) and (12) 
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are « 1, equation (13) can be further 


Since both - — — and - — 
g g 

simplified and the relation between density and temperature can 
be expressed in an approximate form 
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The physical significance of this appro formula is 

that it gives us the order of magnitude of density amplitude 

at the turbopause level once the temperature amplitude is 
specified. However, it does not contain any information on the 
temperature amplitude itself. The temperature amplitude, 
therefore, has to be treated here as the principal parameter 
which, along with the temperature phase, has to be assumed in 
a plausible way. Eqn . (14) furthermore tells us that temperature 
and density at the turbopause should be expected to be in phase. 

In the region between 70 - 120 km, the main source of 

heating is the absorption of the Lyman a and ultra violet 
radiation in the Schumann continuum range and thus the tem- 
perature has a strong diurnal component . Considering only 
the first harmonic terms, the diurnal variations in temperature 

and density may be written in the following form 


T = T o + T i cos (<d t + 0) 


(15) 


p(z Q ) = p o^ Z o^ + MV COS ( fjjt + ^ d 6 ) 

The variations between density and temperature at the turbopause 
according to equation (14) may be related by the following 
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express ion 
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Since (— ) and ^ — y arc' the relative amplitudes for the 

temperature and density, it is evident from equations (14) 

and (17) that for a given temperature variation the density 

variation is amplified by a factor 22. A simple estimate for 

n 

this amplification may bo made as following. If the lower and 
upper boundaries in this problem are assumed to be 70 and 
120 km, respectively, then, assuming a mean scale height of 


zo 

about 10 km, y=— ^ 5. 

T-, D-. , x 

Thus , iX ~ 0 .1 , -0.5. 

The example given above is a rough estimate and is meant 
to be only illustrative. The actual amplification may vary 
over a wide range if the effects of the viscosity and eddy 
diffusion are taken into account . The discussion of this 
section is meant to emphasize the point that a large diurnal 
variation in density at the turbopause level is possible even 
though the temperature variations may not be very large. In 
the next section we shall discuss the effect of this variability 
on the solutions of the heat conduction equation in the diffusive 
region . 
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HEAT CONDUCTION EQUATION 

The numerical results on density and temperature presented 
in this section are based on the solution of the heat conduction 
equation as discussed by Stubbe (1970). (To be referred to as 
ST 1.) For expediency, we have neglected the energy input from 
the electron-ion gas and the convection term arising from the 
horizontal transport. An exclusion of these terms from the 
heat conduction equation may cause the temperature to be under- 
estimated by 10 - 20%. However, this can easily be compensated 
by suitably adjusting the heating efficiency of the incoming 
radiation which is an undetermined parameter in this problem. 
Both the phase and the day to night ratio of the exospheric 
temperature are relatively independent of the choice of the 
heating efficiency, the horizontal convection and the energy 
input from the ionized gas. The results of this paper can, 
therefore, be discussed without any loss of generality. 

With the assumptions outlined above, the heat conduction 
equation for the diffusive region may be expressed in the fol- 
lowing form: 


1 

St z Sz k S c p4 


3 ( . BTv p 

T. ( 1 s¥> h p < 



(18) 


where n^ is the density of the £th constituent of the neutral 
atmosphere and c ^ the corresponding specific heat at constant 
pressure. The neutral constituents of the atmosphere are as- 
sumed to be 0 2 , N 2 , 0, H e and H. The other terms in equation 
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(18) arc* defined as follows: 
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where T (-\) is the exospheric photon flux intensity and nr.( the 
absorption cross section. The factor of proportionality, s , is 
the efficiency of conversion of solar energy into heat energy. 

As in ST-1, e - 0.4 except for the Schumann Runge continuum 
range (1325 \< 1775 8) for which we assumed 6=0.1. In com- 

puting , sec v which accounts for the solar zenith angle 
variation is replaced by an appropriate Chapman function near sun- 
rise and sun-set as discussed in ST-1 . 

P^ =- Heat loss arising from the infrared radiation emitted by 
thermally excited atomic oxygen (Bates 195J ) . 
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"Breathing" velocity of the atmosphere 
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where H(z) is the density weighed mean scale height and p(zo) the 
total pressure at z - zo. 

The first term on the right hand side of equation (22) is 
the additional term arising from the variation in density and 
temperature at the turbopause. 

DIURNAL VARIATIONS IN DENSITY AND TEMPERATURE 

The numerical solutions of the heat conduction equation 
as outlined in the previous section in conjunction with equations 
(1) , (15), and (17) give the diurnal variations in density and 
temperature above the turbopause. The numerical procedures 
adopted here are exactly the same as discussed in ST-1 . The 
neutral constituents of the atmosphere are assumed to be , N^ , 

0, H 0 and H, all in diffusive equilibrium. Their values are, 
however, specified at different levels. The values of H and He 
are specified at 500 km. As in ST-1 

n(H,500) = 1.5 x 10 7 exp { - ^ i ~ 7 ~ & ~ ^ ' 3 cln ^ (23) 

The expression for n(He) at 500 km is adopted from Keating et al. 
(1970) and is given in the following form 

log n (He, 500) = - 14.552 + 13.104 log T (500) -2 . 019 (log T(500)) 2 - 
10 10 10 

0.4 66 
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where and 


are respectively the solar declination angle and 


the geographic latitude both measured in radians. The last 
term in equation (24) is introduced to account for the winter 
helium bulge as discussed by Keating and Prior. 

In the case of 0 ,^, N 2 and 0, these values are specified 
at 120 Km in accordance with (17) 
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120 km is assumed 


to be of the form 


T (120 Km) - 500 [l +0.1 cos (wt 4- 0] (26) 

The numerical solutions presented here correspond to the 
equator (' 0) , equinox conditions ( 1 ~ 0) and for medium 

solar activity. The solar U.V. flux values are assumed to be 
30 percent higher than those given by Hinterreger et al . (1965). 

Figo 1 shows the diurnal variations in density and tempera- 
ture at 300 Km for R - 0.5 and 0 ~ 0. 

The assumption of 0 = 0 implies that the density and 
temperature at the lower boundary are in phase with their 
maximum values at midnoon and minimum values at midnight. 

R = 0.5 implies a factor of 3 variation in density at 120 Km, 
related to the minimum value. The main features of Fig. 1 may 
be summarized in the following: 

1. The density and temperature at 300 Km do not reach 
their maximum values at the same time. Whereas the 
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density maximum occurs at about 1500 hours local 
time in good agreement with the satellite drag results, 
the temperature maximum occurs about two hours later 
as inferred from incoherent back-scatter data. 

2. The temperature minimum occurs about two hours later 
than the density minimum in accordance with both the 
satellite drag and back-scatter results. 

3. The minimum to maximum change in temperature is from 
about 800 to 1200, an increase of 50% from day to 
night. This is in very good agreement with the diurnal 
amplitude in temperature inferred from back-scatter 
measurements. 

Since both R and 0 are assumed parameters in this problem, 
it is important to study the effects of their variations. 

In Fig. 2, the diurnal variations in density and tempera- 
ture at 300 km are shown for varying values of 0. The value 
of R is assumed to be 0.3 in all the cases. We note that the 
diurnal characteristics of the neutral temperature remains 
practically unaffected with the variations in 0. The changes 
in the neutral density, on the other hand, are very significant 
both with regard to their phase and magnitude. Both the density 
and temperature are in the same phase with their maximum values 
at 1700 hours when 0 = -90°, i.e. when the density and tempera- 
ture at the lower boundary have their peak values at 1800 hours. 
As 0 -* o, the density peak moves to an earlier hour. The shift 
in the time of maximum is not in the same proportion as the 
changes in 0. For example, when 0 changes from 90° to o° (a 
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phase change of six hours at the lower boundary) , the density 
peak at 300 km moves only from 1400 hour to 1500 hours. 

Finally, when 0 180° (the maximum density and temperature 

at the lower boundary at midnight), the density maximum at 
300 km occurs late in the evening. 

The change 1 in the diurnal density amplitude with the change 
in 0 is quite noticeable in Fig. 2. The day to night varia- 
tion in density is maximum for 0 -90° and minimum for 0 -■ *90°. 

In Fig. 3 are shown the diurnal variations in density at 
200 and 300 km for varying values of R. The phase angle 0 is 
assumed to be zero and R is assumed to vary from 0.0 to 0.5. 

The shift in the time of maximum with the increasing value of 
R is quite noticeable, particularly at 300 km. Larger day 
to night variations in density are seen with increasing values 
of R. This effect is more noticeable at 200 km where the day 
to night ratio are respectively 1.2, 2, and 3 for R - 0, 

0.3 and 0.5. 

There are practically no change in the times of the maximum 
and minimum and very little change in the magnitude of the 
neutral temperature with the variations in R f^ot shown in Fig. 3. 

The diurnal variations in density in the altitude range 
from 200 to 1400 km are shown in Fig. 4. The values of R and 
0 for this case are taken to 0.5 and 0.0, respectively. The 
important points to note in this figure are the gradual shift 
of the time of maximum to later hours with increasing altitude. 

In the altitude range of 200-500 km, the peak density occurs 
between 1400-1500 hours. The time of the minimum also ad- 
vances gradually to later hours beginning with 1:30 hours at 
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200 km. The diurnal amplitude is also height dependent, 
increasing with altitude up to 600 km. Above 600 km, the 
day to night variation becomes less and less pronounced. 

At about 1000 km there is practically no variation from day 
to night. Above 1200 km, the situation is reversed and the 
day-time density becomes lower than the night-time. This 
change in the diurnal characteristics of the neutral density 
above 600 km can be attributed to the presence of hydrogen 
as an important constituent in this height range. Above 1000 km, 
H becomes a dominant constituent. According to Eqn. (23), 
n(H) varies inversely with temperature because of its escape 
from the gravitational field of the earth. 

At present there is no reliable data for the study of the 
diurnal variations above 1000 km. Most of the satellites in 
this altitude range are of very low eccentricity and hence 
not suitable for this purpose. An observational data on 
diurnal variations in the height-range above 1000 km will prove 
very useful in determining the role of hydrogen in this region 
and the processes controlling its escape from the earth’s 
atmosphere . 

CONCLUDING REMARKS 

We have shown how the present discrepancy in the phase and 
the diurnal amplitude of the thermospheric temperature inferred 
from satellite drag and incoherent back-scatter can be reconciled 
by introducing the concept of dynamic diffusion in the solution 
of the heat conduction equation. We have been able to explain 
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most of tin* observed features of the diurnal eluirae t er ist i- * • 
in density and temperature without introducing a second heat 
source* or involving the concept of horizontal transport. 

Rtshbeih (l‘h)7) has suggested the latter as a mechanism tor 
causing tie phase* lag between the density and temperature 
at thermospheric heights. At present, however, this idea is 
purely qualitative and it is difficult to ascertain its merit 
without carrying out some detailed numerical analysis. 

The concept ol dynamic diffusion may appear somewhat 
artificial, since* the atmosphere cannot be arbitrarily divided 
into two regions With a common boundary at 120 km. The choice 
of 120 km, however, is not very critical and the conclusions 
ol this paper are equally valid it the boundary is chosen at 
lower heights. Ideally, one should treat the entire region 
between 70 and 1500 km as one entity without introducing a 
boundary m between. The problem of getting numerical solutions 
in this ease, however, is formidable. Because of the increased 
aii’ mass below 100 km, there is an additional problem in getting 
converged solutions in an acceptable time. The approach adopted 
in this paper is a reasonable compromise. It is simple and 
easy to adopt for the purpose of generating atmospheric models. 
Thus, in computing the temperature from the density obtained 
from satellite drag measurements, the static diffusion model 
can be easily replaced by a dynamic diffusion model as outlined 
in this pap 
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APPENDIX 

Equation (6) and (7) may be derived as follows: From (2), 

( 3 ) » (4) and (5), the first and zero order equations are ob- 
tained in the following form: 
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From A-l and (5) , the expression for the rth harmonic in 
velocity is given by 
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Substituting (A-6) in (A-2) and assuming the atmosphere to be 
horizontally stratified so that all the horizontal gradients 
may be neglected compared to the vertical, the following 
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differential, equation is obtained 
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In deriving A-7, the altitude variations of T and T, have 

o lr 

been neglected for simplicity. The solution of A-7 is simply 
given by 
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If the sun creates the main driving force, a can be taken to be 
equal to jd ; in this case the expression for and ^ can 

further simplified. Since K is of the order of u> [except when 
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for the vertical velocity as given by (7) can easily be 
obtained from (A-8) and (A-6) . 
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Figure 2 . 
300 km . 
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The diurnal variations in density and temperature at 
The effect of phase change at the lower boundary. 
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Figure 3. The diurnal variations in density at 200 and 300 km. 
The effect of the diurnal amplitude change 
at the lower boundary . 






